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Constructing a universal operation scheme for Majorana qubits remains a central issue for the 
topological quantum computation. We study the Landau-Zener-Stiickelberg interference in a Ma¬ 
jorana qubit and show that this interference can be used to achieve controllable operations. The 
Majorana qubit consists of an rf SQUID with a topological nanowire Josephson junction which hosts 
Majorana bound states. In the SQUID, a magnetic flux pulse can drive the quantum evolution of 
the Majorana qubit. The qubit experiences two Landau-Zener transitions when the amplitude of 
the pulse is tuned around the superconducting flux quanta 2e/ft. The Landau-Zener-Stiickelberg 
interference between the two transitions rotates the Majorana qubit, with the angle controlled by 
the time scale of the pulse. This rotation operation implements a high-speed single-qubit gate on 
the Majorana qubit, which is a necessary ingredient for the topological quantum computation. 

PACS numbers: 03.67.-a, 74.50.+r, 74.90.+n 


I. INTRODUCTION 

Majorana bound states (MBSs) which reside in topo¬ 
logical superconducting systems are drawing much atten¬ 
tion both theoretically and experimentally j71-[T7| . The 
so-called topological superconductivity in this context is 
the spinless p-wave superconductivity, in which a Ma¬ 
jorana number can be defined by the Pffafian of the 
Bogoliubov-de Gennes Hamiltonian Q. MBSs are zero 
energy quasiparticles in these topological superconduc¬ 
tors. They are localized at the ends of one dimen¬ 
sional systems[B[, or the core area of the superconduct¬ 
ing vortices in two dimensional systems [dj. Electrons 
in natural systems always have spin, therefore, MBSs 
are predicted in artificial structures such as the inter¬ 
face of a three dimensional topological insulator and a 
superconductor [5,: 17l| . the spin-orbit coupling nanowire 
in proximity to a superconductor [U, Hj, an d the f er ~ 

romagnetic atom chain on top of a superconductor [l6j . 
These systems have a common character that the spins of 
the electrons near the fermi surface are effectively elimi¬ 
nated by the spin-orbit coupling and the Zeeman energy. 
Then, the effective spinless superconductivity, i.e. the 
topological superconductivity, is achieved through the 
proximity effect. 

MBSs in topological superconductors are interesting 
because of their non-Abelian exchange statistics [H, []j|. 
They provide a realistic platform for investigating a sim¬ 
plest example of non-Abelian particles. Besides, MBSs 
are proposed to be useful in quantum computation^. 
They can construct Majorana qubits, which resist to lo¬ 
cal perturbations and can store quantum information for 
a long time. Furthermore, the Majorana qubits can be 
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rotated by the braiding of MBSsfiH]. Importantly, these 
rotation operations are topologically protected. The ef¬ 
fect of the operation is determined by the topology of the 
braiding, not the detailed path. This topological quan¬ 
tum processing provides a possible scheme for the topo¬ 
logical quantum computation. However, the braiding op¬ 
erations are not sufficient to realize universal quantum 
gates [2]. They must be supplemented by non-topological 
operations on Majorana qubits, which have been pro¬ 
posed bv using quantum dots[10, Il2l. superconducting 
qubits [20l I2lj. or microwave cavities [ 23 . 

One promising approach for operating qubits is the 
Landau-Zener-Stiickelberg (LZS) interference[23j]. The 
LZS interference is a standard quantum phenomena in 
quantum mechanics. It occurs in quantum two-level 
systems in which the two levels exhibit avoided level 
crossings. At the two sides of the crossing point, the 



FIG. 1. (Color online), (a) Schematic setup of the Majorana 
qubit, which consists of an rf SQUID with a tunneling barrier 
junction on the topological nanowire. Four Majorana bound 
states are localized at the ends of the wire and the tunneling 
barrier, (b) The energy spectrum of the Majorana qubit as 
a function of magnetic flux <f>, where the eigenstates are up- 
spin (solid line) and down-spin (dashed line) in the psudo-spin 
representation. 
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physical properties of the eigenstates are exchanged. 
Therefore, non-adiabatic Landau-Zener transitions be¬ 
tween the two levels may happen when the system trans- 
verses the avoided crossing under a time varied control 
parameter [24], [25[. The interference between these transi¬ 
tions is called LZS interference (26j. The LZS interference 
is very useful for coherent quantum operation, since it is 
robust to certain noises a nd p ossible to implement high 
fidelity quantum gates (H, [27(. In Ref. [27], the LZS in¬ 
terference has been experimentally achieved for realistic 
quantum control on quantum dot qubits. 

In this work, we propose to use the LZS interference 
to implement quantum operations on Majorana qubits. 
For this purpose, we adopt the system sketched in Fig. 
1, which consists of an rf SQUID with a topological 
nanowire Josephson junction. The nanowire hosts MBSs 
near the tunneling barrier and the ends of the wire. These 
MBSs form a Majorana qubit, which is correlated with 
the superconducting phase difference across the junction. 
Because of this correlation, the Majorana qubit is de¬ 
scribed by the Josephson energy of the junction, which is 
a 2 x 2 Hamiltonian with avoided level crossings. If a mag¬ 
netic flux pulse is applied, the Josephson Hamiltonian 
will evolve and transverse through these crossings. We 
consider a triangular pulse, for which the maximum flux 
value is tuned around the superconducting flux quanta 
<F 0 = 2 e/h. Within one pulse, the Josephson Hamilto¬ 
nian transverses through the same crossing point twice, 
and then come back to the initial Hamiltonian. In con¬ 
trast, the Majorana qubit state does not come back to 
the initial state after the pulse. Landau-Zener transi¬ 
tions happen at the avoided crossings, and the resulted 
LZS interference rotates the Majorana qubit, where the 
rotating angle is controlled by changing the time length 
of the pulse. Our work provides a one-qubit gate on the 
Majorana qubit, which is important for realizing topo¬ 
logical quantum computation. 

We organize this work as follows. The model for the 
system is presented in section II, then we study the LZS 
interference analytically under this model in section III. 
We present a numerical simulation in section IV. Finally 
we give discussions and a conclusion in section V. 


II. MODEL 

The system illustrated in Fig. la is constructed by 
a nanowire Josephson junction which hosts four MBSs. 
Two of them ( 71 , 74 ) are localized at the two ends of 
the wire, and the other twos ( 72 , 73 ) at the two sides of 
the tunneling barrier. The Majorana qubit built by the 
two MBSs near the tunneling barrier is directly described 
by the Josephson energy of the junction, which is differ¬ 
ent from conventional junctions in two aspects. First, it 
is actually a matrix for a Majorana junction^, rather 
than a pure number for conventional junction. The basis 
states of the matrix are the Majorana qubit states. In 
this sense, it is more accurate to be named as Josephson 


Hamiltonian. Second, the Josephson Hamiltonian has 
the 47 t period components in its diagonal elementsjjl, Q, 
different from the 2-7T period Josephson energy in con¬ 
ventional junctions. With these insights, the Josephson 
Hamiltonian for this Majorana junction has been given 
as @, 

H = *.Ej 7273 cos(7T$/d> 0 ) + *Jl7i72 + *Jr7374, (1) 

where 71,2,3,4 are the four MBSs, Ej is the Josephson 
energy due to the coupling of the two MBSs around the 
tunneling barrier; Sl,r represent the coupling between 
the distant MBSs at the left and the right side of the bar¬ 
rier, respectively; $ is the applied magnetic flux through 
the SQUID, which controls the phase difference across 
the junction in the vanishing inductance regime under 
current consideration. 

We define two fermionic operators fl = (72 + 173 )/2 
and f\ = ( 74 +* 7 i )/2 to construct a fermionic representa¬ 
tion for the MBSs. Then the Hamiltonian is transformed 
to, 

H =Ej( 1 - 2 /}/i) cos ( 7 T$/$ 0 ) + M/2/1 + /I/2) 

+ J_(/ 2 t / 1 t + / 1 / 2 ). 

where = Sl ± Sr. We find that the Josepshon 
Hamiltonian is depending on the occupation states of the 
two fermionic operators, which is nothing but the par¬ 
ity states built by the four MBSs(22j. We are studying 
the quantum coherent evolution of the Majorana states, 
which occurs at a low temperature where the supercon¬ 
ducting energy gap becomes large enough. Thus we 
can ignore the quasiparticle poisoning from high energy 
quasiparticles which strongly suppressed by the super¬ 
conducting gap at low temperature. In this case, the 
total fermionic parity is conserves in this system. The 
Hilbert space of the MBSs can be divided into two dis¬ 
connected subspace with even and odd total fermionic 
parities. Without losing generality, we choose the even 
total parity subspace. In this subspace, we can define 
psudo-spin states |t) = | 0 ) and ||) = /|/|| 0 ), with | 0 ) 
the vacuum states for the two fermionic operators. These 
two psudo-spin states represent the two eigenstates of the 
Majorana qubit. Therefore, a spin representation for the 
Majorana qubit is established. In this representation, the 
Hamiltonian of the Majorana qubit is rewritten as, 

H = £ , jcos( 7 T$/$o)cr 2 +8<J X , (3) 

where a X:Z are Pauli matrices acting on the pusdo-spin 
states, and S = <5_. The two eigenvalues of H are E± = 

±\JEj cos 2 (^) + S 2 , which are depicted in Fig. lb as 
a function of magnetic flux <f>. This Hamiltonian is a typ¬ 
ical two-level qubit system with avoided level crossings. 
Away from the crossings, the coupling between spin-up 
and spin-down state is small compared to the energy dif¬ 
ference. Then these two spin states are the eigenstates 
of the qubit. The qubit dynamics is adiabatic. How¬ 
ever, the energy difference between the two spin states 
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vanishes at the crossing points, then the spin coupling 
5 dominates the dynamics of the qubit states. When 
the system traverses through the crossings under a time 
varying magnetic flux, interesting physics occurs: non- 
adiabatic Landau-Zener transitions may happen, and the 
qubit state will experience a rotation on the Bloch sphere 
representation( 23 j. This Landau-Zener transition has 
been extensively studied in the literaturc[ 23 j, with real¬ 
istic examples achieved in atomic systems, quantum dot 
systems, and superconducting systems. 


III. LANDAU-ZENER-STUCKELBERG 
INTERFERENCE 

The Landau-Zener transition can significantly modu¬ 
late the quantum state of the Majorana qubit. However, 
the system will not restore to its original Hamiltonian af¬ 
ter one Landau-Zener transition. The Hamiltonian will 
be added with magnetic flux, which brings an obstacle 
for the realistic quantum control on the Majorana qubit. 
One solution to this obstacle is to consider a magnetic 
flux pulse. During the pulse, the flux first increases and 
the system transverses through an avoided crossing point, 
then the flux decreases and the system transverses back 
through the same crossing point. Finally, the flux van¬ 
ishes after the pulse and the system will restore to the 
original Hamiltonian. The Majorana qubit transverses 
through the same crossing point twice within one pulse. 
Two Landau-Zener transition happens and the LZS in¬ 
terference will rotate the Majorana qubit. 

A flux pulse with a triangular shape is a natural choice 
for this purpose. It is one of the simplest pulse shapes, 
and provides a constant increasing and decreasing speed 
for the magnetic flux. This linear dependence simplifies 
the analytic solutions. We consider an triangular pulse 
with a function of, 
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FIG. 2. (Color online), (a) The asymmetric triangular flux 
pulse as a function of time, (b) The time evolution of the 
two diagonal terms of the Hamiltonian Eq. ©. (c) The time 
evolution of the Majorana qubit sate represented by the two 
components |i/’i(t)| 2 (solid line) and \tp 2 {t)\ 2 (dashed line). 
The parameters are A = 1, 8/Ej = 0.001, toi/8 = 0.2 and 
102/8 = 0.025. 


crossing $ = $0/2. Near t\ and £2, the diagonal par t of 
the Hamiltonian is small and can be linearized as[ 23 j, 


e{t) = Ej cos(tt$(/)/$ 0 ) ~ ±Vi ,2 (t - * 1 , 2 ), (5) 


where 


Vl,2 = Ej 


d_ 

dt 


cos(7T$(/)/$ 0 ) 


t=t 1,2 




( 6 ) 


Then the effective Hamiltonian is linearized near the 
crossing times, 


|W, 0 <i<^ 

<f>(/) = $ 0 x \ A-w 2 (t- £), £<t<£ + £ (4) 

I 0, others 

where the amplitude A and the velocity to 1.2 are posi¬ 
tive. As shown in Fig. 2 a, this pulse starts from time 
zero, increases with a constant speed to 1, reaches a max¬ 
imum value of A$o, then decrease with a constant speed 
of L02 , finally vanishes at time A/to\ + A/to 2. This tri¬ 
angular pulse is asymmetric for general case. It becomes 
symmetric under the special condition of ui\ = to 2. In the 
following calculations, we find that the asymmetric pulse 
is more suitable for our purpose to control the Majorana 
qubit. 

Landau-Zener transitions appear in this system only 
when the pulse amplitude A is greater than 1 / 2 . For sim¬ 
plicity, we consider a pulse amplitude ofl/ 2 <A< 3 / 2 . 
Then only one avoided crossing point is involved dur¬ 
ing the pulse. In the following analysis, we denote 
ti = l/( 2 wi) and <2 = t\ + (A — 1/2 )/w2 as the cross¬ 
ing time when the system transverses across the avoided 


H lt 2{t) = Vi t 2(t-ti l 2)a z +Sa x . (7) 


These linearized Hamiltonian H\ 2 are typical examples 
of Landau-Zener problem. The Landau-Zener transition 
probability for a single crossing is given bylMiH 

-Pi,2 = exp(- 27 r/ 3 i j2 ), (8) 

where 61,2 = d 2 /2v\p. If the system starts from the 
lower level initially, then the Landau-Zener transition 
probability P± t 2 describes the probability of upper level 
occupation after one crossing. As shown in Ref. [ 28 ], 
these Landau-Zenner transitions can be described ap¬ 
proximately by an unitary evolution matrix Ni 9 [ 23 ll 28 j. 


Nl,2 


Jl\2 V 1 ~ \ 

( 9 ) 


where tp i j2 = -f + /3i,2(ln/3i j2 - 1) + argT(l - */?i >2 ), T 
is the gamma function and the phase tpi.2 is monotonous 
function changes from 0 in the adiabatic limit (Pi j2 —>• 
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0 ) to 7r/4 in the diabatic limit (Pi .2 —> These 

two evolution matrices connect the wave function of the 
system before and after the Landau-Zener transitions, 
which reads \I/(ii )2 +0) = Ni, 2 ^{ti ,2 — 0) with 'h = |t 

} + ^2 14-)- 

Now we consider the LZS interference between the two 
Landau-Zener transitions, and obtain the finial Majorana 
qubit state. The system evolution between these two 
Landau-Zener transitions is approximated by an adia¬ 
batic evolution matrix [2^, (28|, 




e -iC(pi) 0 \ 

0 e iC(t '’ 4) ) ’ 


( 10 ) 


where £(/' ,t) = ^ J* cos ^p-dr are adiabatic evolu¬ 
tion phase. The quantum state rotation induced by the 
LZS interference is written as [13, 

*(i) « U(t, t 2 )N 2 U(t 2 , t^NiUttuto^o, (11) 

where di'o = (0,1) T is the initial state for the Majorana 
qubit, and ’L(i) = 2 (t)] J is the final state. We 

take U(ti,to) as identity matrix with appropriate choice 
of the phase for the basis states. Then the final qubit 
state reads, 

ipi(t)= [\M(1 -P 2 )e if(t2 ’ tl)+iv ~ 2 (12) 

+V(1 - P 1 )P 2 e- iat2 ’ tl)+i ^} e - iat ’ t2) . 

We immediately find that the Majorana qubit state has 
been rotated. The phase of ipi(t) experiences the Larmor 
precession, therefore will rotate after the pulse horizon¬ 
tally in the Bloch sphere. In principle, any horizontal 
rotation angle for the Majorana qubit can be achieved 
with appropriate precession time[27[. Here, we concen¬ 
trate on the longitudinal rotation of the Majorana qubit, 
which is given by the amplitude of the final wave func¬ 
tion, 


|Vh| 2 = (Pi+P 2 - 2 P!P 2 ) ( 13 ) 

+2\/PiP 2 (l - Pi)(l - P 2 ) cosy, 

where % = 2 ((t 2 ,ti) — (p\ + <^ 2 comes from the adia¬ 
batic phase evolution, which is sensitive to the parame¬ 
ters of the pulse [HI- The first term in Eq. m provides 
the rotation angle of the Majorana qubits from (0,1) to 
[ipi(t),ip 2 (t)\. The rotation angle is fully determined by 
the probability of the two Landau-Zener transitions Pi 
and P 2 , therefore is controlled by the flux varying speed 
ui\ and u> 2 . The second term can be treated as the uncer¬ 
tainty for the rotation operation, since the phase factors 
are uncontrollable. We could reduce this uncertainty by 
modulating Pi or P 2 , making one of them approaches 0 
or 1. This regime can be achieved only for the asymmet¬ 
ric triangular pulse, since the symmetric triangular pulse 
leads to Pi = P 2 . Then the uncertainty term will always 
be significant. 
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FIG. 3. (Color online) The time evolution of the Majorana 
qubit sate represented by the two components \ij)\{t)\ 2 (solid 
line) and \i/) 2 (t)\ 2 (dashed line), with (a) L 02 /5 = 0.025, (b) 
UJ 2 /S = 0.01, (c) uj 2 /S = 0.004, and (d) lo 2 /S = 0.002. Other 
parameters are A = 1, S/Ej = 0.001, and uii/5 = 0.5 
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IV. NUMERICAL SIMULATION 


Now we study the system described by Eq. (3) nu¬ 
merically to obtain quantitative results. First we setup 
the equations for numerical simulations. We expand the 
Hilbert space of the Majorana qubit system with the 
psudo-spin states |\H(t)) = |f) + 1 ^ 2 (t) |4-)- Then 

the Schrodinger equation for the Majorana qubit is writ¬ 
ten explicitly as, 


i— ( ^ 

dt V ^2 


E J cos ^ 

6 -Ej cos J V ^2 ) ' 


(14) 

Numerical simulations for this Schrodinger equation are 
performed with standard finite difference method, in 
which the evolution operation is linearized within each 
small segment of time. We consider the asymmetric tri¬ 
angular magnetic flux pulse described in Eq. (4). The 
signal starts from t = 0 and ends at t = A/oji + A/oj 2 , as 
shown in Fig. 2a. This pulse increases the flux through 
the SQUID from zero to A$ 0 rapidly, then reduces back 
to zero slowly. During this process, the two diagonal 
terms of the Hamiltonian, which is shown in Fig. 2b, 
variate and cross with each other twice. Around these 
two special crossing points, the diagonal part vanishes, 
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thus the off-diagonal part dominate the system and the 
qubit experiences Landau-Zener transitions. We calcu¬ 
late the wave function evolution of the Majorana qubit 
starting from = (0,1) T , and show the results in Fig. 
2c. The qubit is staying at the initial state before the 
pulse. After t = 0, the pulse is applied and the Hamilto¬ 
nian begins evolution. We find that the system stays at 
the initial state away from the crossing points, because 
the energy difference between the diagonal part of the 
Hamiltonian is big enough to induce adiabatic dynamics. 
When the flux increases and the system approaches the 
avoided crossing, the evolution becomes non-adiabatic, 
and the qubit will have a probability of psudo-spin ro¬ 
tation. This is the Landau-Zener transition as expected 
in theoretical analysis. The increasing speed of the flux 
is very quick for this pulse. Therefore, the evolution is 
in the diabatic limit. The psudo-spin rotation is small. 
When the flux decreases, the second Landau-Zener tran¬ 
sition happens. This time, the decreasing speed of the 
flux is moderate. The LZS interference between these 
two transitions then rotates the qubits, which is clearly 
demonstrated with the two components of the wave func¬ 
tion in Fig. 2c. 

This qubit rotation is controllable by modulating the 
time scale of the magnetic pulse. More specifically, it is 
controlled by the flux decreasing speed of the pulse, since 
we increase the flux fast enough to achieve the extreme 
diabatic behavior. In this scenario, the error term due 
to the phase accumulation in the adiabatic region is re¬ 
duced. We show the results with different rotation angles 
for the Majorana qubit in Fig. 3. The rotation angle is 
determined by the flux decreasing speed of the pulse ui 2 - 
We find that the rotation can be as large as an inversion 
from (0,1) state to (1, 0) state. With these rotation op¬ 
erations, we achieve a one-qubit gate for the Majorana 
qubit simply by applying a magnetic flux pulse. This 
one-qubit gate is a good supplement to the braiding op¬ 
erations, and should be important for realizing universal 
quantum gates in topological quantum computation. 


V. DISCUSSIONS AND CONCLUSION 

In this work, we exploit the asymmetric triangular 
pulse to implement the one-qubit control. The trian¬ 
gular pulse has the advantage of being simple in theory, 
yet is experimentally easy to achieve. However, a sym¬ 
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